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Abstract. In this paper we consider flow-equations where we allow a normal ordering 
which is adjusted to the one-particle energy of the Hamiltonian. We show that this 
flow converges nearly always to the stable phase. Starting out from the symmetric 
Hamiltonian and symmetry-broken normal ordering nearly always yields symmetry 
breaking below the critical temperature. 
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1. Introduction 

In applying flow equations to a Hamiltonian |T] one typically starts out from a 
Hamiltonian which does not show an explicit symmetry breaking even if the symmetry 
of the system will be broken below some temperature. As in the calculations on the 
Hubbard model [21 El HI E] one can direct the flow to a form in which molecular-field 
approximation becomes exact. This effective Hamiltonian is still not symmetry broken. 
Only the molecular-field formalism breaks the symmetry below the critical temperature 

T 
- 1 - ci- 
In the fermionic renormalization group flow [HI El 13 ED] the vertex functions will 
at least within weak-coupling approximations diverge at some length scale below T c , so 
that the approximations become unreliable and one has to resort to other methods in 
this regime. Thus it is desirable to have a way to introduce symmetry breaking from the 
beginning. Recently Salmhofer, Honerkamp, Metzner, and Lauscher ^2] have added a 
symmetry breaking field to the Hamiltonian and showed that this additional field leads 
into the symmetry broken phase. 

For the Hamiltonian flow it is not necessary to add a symmetry breaking term to 
the Hamiltonian. Instead it is sufficient to choose a normal-ordering which is symmetry 
broken (compare JI]). One can show that the system will nearly always converge to the 
stable state, that is in case of symmetry breaking (that is below the critical temperature) 
it runs to a symmetry broken state, whereas if the symmetric state is stable (above T c ) 
then it will run to the symmetric state. 

With this approach the normal ordered form of the Hamiltonian shows an explicit 
dependence on the normal ordering. We adjust continuously the normal ordering to 
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the one-particle contribution of the Hamiltonian and we can show that the normal 
ordering has the interesting property to run in the direction of the stable state. This 
implies that starting out with a symmetry-broken normal ordering, then above the 
critical temperature it will converge to the symmetric state whereas below criticality it 
converges to a symmetry broken state. 

We start out from the flow equation [I] 

where r] is chosen so that the Hamiltonian approaches a diagonal or block-diagonal form. 
We write H in the normal-ordered form 

H =: Hq : g (2) 

where G defines the normal ordering 

Gkj = (afcdj)o- (3) 

Here we denote both creation and annihilation operators by a. Then the flow equation 
reads 

dH G 5H G dG 

^T = ^ h]g + 1gW (4) 

The second term indicates the change of Hq due to the change of G. G itself can 
be described as the expectation values of corresponding to a bilinear Hamiltonian H°. 
Then the flow equation is supplemented by a second equation for H°. This equation 
describes the adaption of H° to the one-particle energy of H, and will be given later 
(HU). 

In the next section we consider the effect of normal ordering on operators. Then we 
give the general relation between a bilinear Hamiltonian and the expectation values in 
thermal equilibrium. Then we return to the flow equation and show that almost always 
the normal ordering will approach a stable state. 

2. Normal Ordering 

The idea behind normal ordering is to subtract expectation values 

G k j = (a k aj) (5) 

from products of operator pairs a k , aj. More precisely one defines normal ordering which 
is indicated by two colons by 

: 1 : =1, (6) 

: ctA(a) + (3B(a) : = a : A(a) : +(3 : B(a) :, (7) 

d 

a k : A(a) : = : a k A(a) : + £ G kj : — A(a) :, (8) 
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where a, (3 are c-numbers. These equations hold for bosons and fermions. In performing 
the derivatives for fermions one has to consider that the operators a anticommute. 
Iteration of (JEJ) yields 

d d 
a kl a k2 ...a km =: (a kl + ^2G kljl -^—)(a k2 + ^2G k2j2 -^—)...a km : (9) 



which can also be written 
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a kl a k2 ...a km =: exp \^G kj right left a kl a k2 ...a km : (10) 
\kj ctafc / 

This is Wick's first theorem The superscripts left and nght indicate that we always 
pick a pair of factors a and perform the derivative on the left factor and the derivative 
t^- on the right factor, so that the factor G k j depends on the sequence of the operators. 
The exponential appears in the equation for the following reason. If we perform the 
operation G-£^ on m pairs of factors a, then there are due to the permutation symmetry 
m\ contributions. Therefore in order to obtain the contribution with factor one we have 
to divide the m-th power of G-£^ by m\, which yields the exponential. Note that for 
fermions the operators a as well as the derivatives anticommute. 

Since in the following we will change the normal-ordering it is appropriate to 
indicate to which expectation values G it is performed. Then we obtain quite general 
for operators A 
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:Ao: =A, A G = exp (E^^^J A(a) (1!) 

If we dissect G k j = Q k j + C k j, where C k] - = Cj k for bosons and C k j = —Cj k for fermions, 
then we have 

: A : G =: A Q : Q , A a = exp \l £ A Q . (12) 

and we need no longer distinguish the sequence of the factors. Suppose we have 

H = v {0) + ^Y. v kj a kaj + 4| H v$ mn a k aja m a n (13) 

kj ' kjmn 

and Vf?j and v k 2 j mn are completely symmetric for bosons and antisymmetric for fermions, 
resp., then the normal-ordered form reads 

u (o) , 1 (i) , 1 (2) (, A \ 

h g =vh + 2i v G,kj a k a j + -^ v G'k jm n a kaja m a nj (14) 

$ = v (0) + \Y, v$G kj + ± Y, v [ $ mn G k p mn , (15) 

kj kjmn 

1 



W _ J 1 ) _i_ 1 ^J 2 ) 



V b[kj — V kj + r, ^L V kjmnGmni (16) 
* mn 

v {2) -v (2) (17) 

U G,kjmn ~ U kjmn- \ L ' I 
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Expressing A and B by A G and B G by means of eq. QlUjl and then transforming back 
to the normal-ordering we obtain for the product of two operators 

: A G (a) : G : B G {a) : G =: exp G kj -^-\ A G {a)B G {b) : \ G<b=a . (18) 
3. Expectation values for H° 

In the following we will consider only fermions. A conventional way to introduce normal 
ordering is to use a Hamiltonian 

H° = ^ikjdk'Cij (19) 
where we use the notation 

4 = a k *, a k = a\*. (20) 
Considering a a column vector and a) a row vector we may write 

e hJ a k .a j = ( o) a T ) I ( ° ) = ( o> a T ) ( * _* r ) ( ° ) 

A f = A, B T = -B. (21) 

Thus e has the properties 

et = c - r gr = - = ( J J ) (22) 

e can be diagonalized with diagonal matrix elements e k and — e k . (This diagonalization 
is performed by a canonical transformation which is isomorphic to a real orthogonal 
transformation. This can be easily seen if one introduces the hermitean linear 

combinations q k = and p k = ^ a */^ .) 

In thermal equilibrium one obtains for diagonal e, i.e. e k j = 5 k j€ k 

so that even when e is not diagonal 

G = si — = 1-1 tanh A (24) 
e^ e + 1 2 2 v 2 ; K J 

holds with Gfcj = G k *j. A variation of e yields 

- — / dre^ee^ 6 "-^ — . (25) 
; + lio e^ e + l v ; 

Thus we may write 

3G k j = —T k j tPq Se P q (26) 

with 

r ^=£Htt%(^?hi)< (27) 



SG = — 8eP* -== = - , 
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4. Free Energy and Stability 

It is well known, that for a given Hamiltonian H and temperature T the free energy 
assumes its minimum for the corresponding statistical operator p = e~^ u jZ. Thus 
one often determines approximately the free energy for the statistical operator p° = 
e -P H ° jz° with H° bilinear in the operators a and determines e so that the corresponding 
free energy F° becomes minimal. One obtains 

F° = E-TS, E = v { q\ S = -ytr(GlnG + (1 - G) ln(l - G)) (28) 

We use that (A(a))° = AG\ a=0 - We will vary this expression. By means of 

^ = A + \^5 GliI ^-A a (29) 

+ i Y. iG kl 6G m f A G + 0{SG'-) (30) 
o kjmn oajOa k oa n oa m 



we obtain 

1 

2 

A variation of S yields in first order in SG 



E = 4 0) + ~A\^G k3 + -v%^ n 8G kj 8G mn + ... (31) 



5S=-^\n(-^0\ SG kj (32) 



jk 

Therefore we obtain in first order in 5G 



SF° = \ \v%\, 3 + k B T\n (^Jj ^ ] 6G kj 



-(v {1) 



e kj )8G kj . (33) 



In order that F° is an extremum one has to choose e k j = v Gk *j, that is the one-particle 
contribution of our Hamiltonian Hq has to agree with H°. In second order in SG we 
obtain 

SF° = lvg* jm *JG kj 5G mn - ~6~e kj SG kj , (34) 
where 5e is given in terms of 5G from (|2*Bj) 

fi^kj = (-T )kj,mnfiG mn . (35) 

which yields 

5F° = ^4, jm , n + (T-% hmn )5G k3 5G mn . (36) 

Only if this expression is positive definite then it corresponds to a stable solution. We 
note that V is positive definite and therefore T _1 exists and is positive definite, too. 
This can be seen if we switch to the basis in which e is diagonal, e k j = 5 k j€ k . Then 

^kj,pq = $kpfijq£ \ji (37) 
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Thus T is diagonal with positive matrix elements along the diagonal. 

These expressions are closely connected to the response function in random phase 
approximation. Suppose one adds a perturbation £t>( lext ) to the Hamiltonian H which 
produces a change 5G ext in the expectation values (a^j) then we have from eqs. (|33I36|) 

\ + (\4\* jm *n + (r-%, m „)5G mn ) SGff = 0. (39) 

If we now relate 5G cxt to a perturbation <W 1,eff ) which would have the same effect without 
the interaction according to eq. (|26|) 

3G mn = —^mn.pq $ V pq ^ (40) 

then we obtain 

r (l.ext) _ ,1 (2) -p . r (l.eff) 

0U G,k*j — ^2 U G,k*jm*n L mn,pq + Ofcp Ojq) OV p , g l^lj 

This factor ^v^T + 1 enters for example in the Lindhard expression for the static 
dielectric constant. One has only to introduce v^ 2 ' oc e 2 /q 2 and 

r d _ f( e i) ~ f( e k) (A0 . 
kj ~ e k - e 3 (42j 

with the Fermi function /(e). The dielectric constant has to be positive for a stable 
system. 

5. Flow Equations 

Now we have to introduce our flow equations. We have two contributions to the change 
of He-, one from the generator of the flow, 

[r 1 ,H] G = g^+ l -g^a k a ] + ... (43) 

and one from the change of the normal-ordering which yields 

ldG k j d 2 1 de mn ( (i ) 1 (2 ) \ , , 

Thus we obtain the change of 

d v Vj _ (i) _ l r (2) de mn 
qj — 9k* j 2 Pii mnV v*qk*j v 40 / 

for the one-particle energy. On the other hand we wish to adapt e so that it approaches 
v^ 1 '. Therefore we introduce a flow equation for e 

= 7(«£} " hi) (46) 

with some positive constant 7. Evidently for fixed the energy e approaches 
exponentially From both equations we obtain 

d( V k*j ~ hj) _ (1) _ /r r Ip (2) w (1) _ ~ s (A7 x 

9k*j lyOkmOjn T 2 PI > mn ^ P* qk* j ) \^ m* n ^mri) V 4 '/ 
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Thus as g>W from the generator of the flow decreases also — e will decrease provided 
the kernel 1 + ^Tv ^ is positive definite. Otherwise the difference i>W — e will nearly 
always increase exponentially. The condition that this kernel is positive definite is 
equivalent to the stability obtained from F° in eq. (|3b'|) above, since the kernel for F° 
differs only by the factor T, which itself is positive definite. 
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